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Abstract 

We show that the path integral for the three-dimensional SU(2) BF theory 
with a Wilson loop or a spin network function inserted can be understood 
as the Rovelli-Smolin loop transform of a wavefunction in the Ashtekar con- 
nection representation, where the wavefunction satisfies the constraints of 
quantum general relativity with zero cosmo logical constant. This wavefunc- 
tion is given as a product of the delta functions of the SU(2) field strength and 
therefore it can be naturally associated to a fiat connection spacetime. The 
loop transform can be defined rigorously via the quantum SU(2) group, as a 
spin foam state sum model, so that one obtains invariants of spin networks 
embedded in a three-manifold. These invariants define a fiat connection 
vacuum state in the q-deformed spin network basis. We then propose a mod- 
ification of this construction in order to obtain a vacuum state corresponding 
to the fiat metric spacetime. 
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1 Introduction 



In a recent review of loop quantum gravity [1], Smolin has given a compre- 
hensive description of a vacuum state for a quantum general relativity (GR) 
theory with a positive cosmological constant A. This state can be defined 
in the loop representation as a linear combination of the q-deformed spin 
network states, where the coefficients are given by the loop transforms of 
the Kodama wavefunction. These coefficients can be understood as the in- 
variants of embedded spin networks in the spatial manifold E determined by 
the path integral of the SU (2) Chern-Simons (CS) theory. These invariants 
arc generalizations of the Jones polynomial [2], and can be expressed as the 
Kaufman bracket invariants, i.e. as evaluations of the q-deformed spin net- 
works [3, 4], or more generally as the Reshetikhin-Turaev (RT) invariants 
[5]. 

However, there exists another set of spin network invariants for the three- 
dimensional manifolds, which are also determined by a path integral over 
the connections, but instead of the CS action one uses the BF theory action 
[6, 7, 8]. The purpose of this paper is to demonstrate that these new invari- 
ants can be understood as the loop transforms of a flat connection vacuum 
wavefunction which corresponds to quantum GR with zero cosmological con- 
stant. This vacuum wavefunction is given by a product of the delta functions 
of F, where F is an SU (2) two-form field strength, and therefore that wave- 
function is annihilated by the GR constraints in the Ashtekar connection 
representation. One can then naturally associate this wavefunction to a clas- 
sical vacuum solution with F — There are many such solutions, and we 
will be interested in constructing a quantum state describing the flat space- 
time metric solution. Our approach is an extension of a partially reduced 
phase space quantization approach of Baez [9], where he imposes classically 
the F = constraint and then considers the wavefunctions on the moduli 
space of flat connections. 

In section 2 we show how the BF theory path integral for a spin network 
can be represented as the Rovelli-Smolin loop transform of a flat connection 
vacuum wavefunction in the Ashtekar connection representation. We then 
show how the constraints of quantum GR for zero cosmological constant 
annihilate this vacuum wavefunction, and how to construct the corresponding 
vacuum state in the spin network basis. In section 3 we define the BF theory 
spin network path integrals as the spin foam model state sum invariants 
for the quantum SU (2) group. In section 4 we present our conclusions and 
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discuss our results, including a discussion on the relationship between the 
F = classical solutions and the flat spacetime metric solution. We then 
propose a modification of our construction in order to obtain a flat spacetime 
metric vacuum state. 



2 Ashtekar formulation 

The Ashtekar formulation of GR [10], can be understood as a special choice 
of variables in the canonical formulation of GR based on the triad one forms 
e"" = e'^dx'^, where a = 1,2, 3, and are the coordinates on the spatial three- 
manifold S. The three-metric is given by hij = e°;eja and the spacetime 
manifold is given by E x R. Let be the inverse triad components, and let 
pf be the corresponding canonically conjugate momenta. Then the Ashtekar 
variables can be defined as 

At = u;t{e)+z^ , F^ = v^e^ , (1) 

where the coefficient i — \/— 1, h — det(/ijj) and the a;(e) is the torsion-free 
SO{3) spin connection. Hence the new connection A is complex, and the 
advantage of this unusual choice is that the Hamiltonian constraint becomes 
polynomial in the new canonical variables {—iA, E), i.e. 

Ch = €abc {E^'E^'F^j + A€ijkE"'E^''E'"') , (2) 

where F is the curvature two form associated to the SU (2) connection A 
and A is the four- dimensional cosmological constant. The Gauss and the 
diffeomorphism constraints are also polynomial in the new variables 

Ca = d,El + ea,,A^E^^ , Q = F^.Ei . (3) 

Because of this simplification, one can adopt a Dirac quantization strategy 
where the pair {A, E) is represented by the operators 

= A" , Ei^-— , (4) 

acting on the wavefunction(al) ^{A), such that 

^^ = , Q*-0 , Ch^^^ , (5) 
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where the C operators are obtained by substituting (4) in the classical ex- 
pressions in a chosen order. The SU{2) connection A is complex, and the 
reality properties of the expressions (1) require that the wavefunction '^{A) 
should be a holomorphic functional of A, i.e. independent of A [11]. This 
holomorphy also implies that in order to construct a wavefunction of the 
triads e, it is sufficient to consider only the real A''s [11]. 

Hence by constructing a gauge invariant and a diffeomorphism invariant 
functional of the SU{2) connection A defined on the three-manifold S, one 
will solve the Gauss and the diffeomorphism constraints. Therefore the prob- 
lem is reduced to analyzing the action of the Hamiltonian constraint operator 
Ch- For example, the Kodama wavefunction 

^k[A] = exp ( TriAAdA+'^AAAA A)] = e^^^^I^l , (6) 
V2A7s 3 / 

which is an exponent of the Chern-Simons action, is gauge and diffeomor- 
phism invariant functional and satisfies the Hamiltonian constraint for a non- 
zero cosmological constant [12]. 

A more advanced strategy of quantization is to use the Rovelli-Smolin 
loop representation [13], i.e. if we think of '^{A) as then one can intro- 

duce the loop states I7) associated to the Wilson loop functionals [7, A] = 
Tr P exp A, and define the loop representation as 

^(7) = = j VA{^\A){A\^) = j VAW[^,A] ^{A) . (7) 

If ^(A) is a solution of the gauge and the diffeomorphism constraints (3), 
then ^^(7) will be a three-manifold invariant of the loop 7. More generally, 
one can replace the 7 with a labelled graph F, and hence obtain an invariant 
of the embedded spin network. In the case of the Kodama wavefunction, 
one can argue that ^(F) should be the invariant given by the Kauffman 
bracket [14] ^ This is a consequence of the fact that the CS action defines a 
topological theory on S, and the path integral (7) defines an observable of 
the theory [2] . 

^This argument is straightforward in the case of the EucUdian GR, when the Ashtekar 
connection is real and the Kodama wavefunction is exp(zA:S'c5). In the Lorentzian case, A 
is complex, and the Kodama wavcfmiction is exp(fcS'cs), so that a factor of i is missing. 
There is no rescaling A cA, where c is a complex number, which will give that factor, so 
that the argument is that a Wick rotation fc — > ifc of the Euclidian result has to be done. 
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However, there exists another topological theory in three dimensions, 
closely related to the CS theory, namely the BF theory (for a review and 
references see [15]). It is also a theory of flat connections, but its phase space 
is bigger, because it includes a one form B field. The BF SU (2) theory can be 
considered as the three-dimensional Euclidian GR with a zero cosmological 
constant [15]. The partition function of this theory 

= JvAVBe'I , (8) 

can be defined via simplical decomposition of the manifold as a spin foam 
state sum [16, 17, 15]. This gives the Ponzano-Regge state sum model[18]. A 
topologically invariant way to regularize the Ponzano-Regge state sum is via 
the SU{2) quantum group, so that one obtains the Turaev-Viro state-sum 
model [19]. 

Given these properties, one is led to consider the following quantity 

=|p^PSiy[7,A]e^/^'-(^^^) , (9) 

i.e. the BF theory analog of the CS Wilson loop expectation value. This 
quantity should be gauge and diffeomorphism invariant, and it can be defined 
via simplical decomposition of E as a spin foam state sum [6, 7, 8]. This state 
sum can be then regularized in a topologically invariant way via the quantum 
SU{2) group. 

The form of the path-integral ^'o is such that after performing the B 
integration we obtain 

^^0(7) = jvAW[-i,A]5{F) , (10) 

where 

5{F) = / VBe'I^'^^^''^ = n Y[6{e'^^^^''^) . (11) 

These are heuristic expressions, which have to be properly defined (see the 
next section), but for the moment if we think of S{F) as a functional of 
the connection A, then S{F) = \E'o(^) is gauge invariant due to group delta 
function identity S{ge^g~^) = S{e^). Furthermore we have 

C^^o-J^Ft^{x)l[me'-''')-0 , (12) 
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and 



.abc 



s s 




(13) 







6Af{x) 54 (x) 
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due to F(x)5(e^(^)) = 0. Hence ^o(^) 

is gauge and diffeomorphism invariant 
functional which satisfies the Hamiltonian constraint for the zero cosmologi- 
cal constant. 

This heuristics imphes that the invariant (9) is a loop transform of a 
quantum gravity F = vacuum state ^o(^) = ^F). This vacuum state 
can be naturally associated to the space of F = classical solutions of the 
Einstein equations. A particular F = solution is a flat metric space when 
A = u{e). Hence understanding the properties of the S{F) state will be a 
first step toward constructing a fiat metric vacuum state. 

One can construct a fiat connection vacuum state in the diffeomorphism 
invariant basis of spin network states |r) [20], as 



(r|0) = *o(r)= VAW[r,A]S{F) ^ VAVBW[r,A]e'I^'^^''^^ 



is a spin network invariant to be defined, and [F, A] is the spin network 
generahzation of the Wilson loop functional, i.e. a product of the edge 
holonomies contracted by the vertex intertwiners. 

3 Spin-foam state sum representation 

Let us now write the spin foam state sum representation of the spin network 
invariant (15) in order to better understand its properties. This path inte- 
gral can be defined via the discretization procedure provided by a simplical 
decomposition of E. Let C be a simplical complex associated to E, and let 
C* be the dual complex. Wc denote by r, A and e a tetrahedron, a triangle 
and an edge of C respectively, and hj v, I and / the corresponding dual cells, 
i.e. a vertex, a link and a face. Let the B one-form be piece-wise constant in 
the 3-polytopes formed by the pairs (e, /), where / = e*, such that the only 



|o) = Eir)(r|o) 



(14) 



r 



where 



(15) 
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non-zero spatial components of B are those parallel to the vectors Then 

/ Tr{BAF) = Y,Tr{B,Ff) , (16) 

where 

Be^jB , Ff = jF . (17) 
Given the variables (17), we can define the path integral (8) as 

z= /n^^^n^^^e^^.^'-^^^^/) , (18) 

I e 

where Ai — J^A. We will fix the measures of integration dA and dB such 
that 

lT{d9i\{5[\{gA , (19) 

I f \im J 

where gi = e^K Note that (19) is a discretized form of the path integral 
/ VA6{F) and therefore JJfd (jli^of 9i) =11/^ {.^^^^ considered as a 

discretized form of the wavef unction S{F). 
By using the group theory formulas 

oo 

S{g) = dim J xAq) = E (2i + 1) xM , (20) 

j 2j=0 

and 

L(2) (5)^^3(5) = (C^igla)* > (21) 

where D^^^ is the representation matrix in the representation j, Xj is the 
corresponding trace and C's are the Clcbsh-Gordan coefficients, or the 3j 
intertwiners, as well as the associated graphical calculus [17, 15], one arrives 
at 

^ = En dimj/n^tet(j/i(r,),.-,J/6M) ■ (22) 

3f f " 

Atet is the tetrahedron spin network amplitude, given by a contraction of four 
3j tensors^, which is also known as the (6j) simbol. The terms in the sum (22) 

■''This is an approach to tlic discretization introduced in [23]. It is a more convenient 
approach for treating non-linear actions then the usual approach where the B field is taken 
to be a distribution concentrated along the edges e [17, 15]. 

§For a tensorial approach to spin networks see [21]. 
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can be interpreted as the amplitudes for a colored two-complex associated to 
C*, which is called the spin-foam [15]. The Ponzano-Regge state sum formula 
can be obtained by by relabelling the sum (22) as a sum over the irreducible 
representations (irreps) which are associated to the simplices dual to the spin 
foam, i.e. the edges and the tetrahedrons of the manifold triangulation 

Z = ^n <iimjenAet(^(r),---,ies(r))=En(2je+l)n(6j)W • (23) 

The expression (23) is not finite, because the infinite sum over the SU (2) 
irreps j divereges. Formally, the sum (23) is topologically invariant, because 
it is invariant under the three-dimensional Pachner moves. It is possible to 
regularize this sum while preserving the topological invariance, by replacing 
the SU{2) irreps with the irreps of the quantum group SUq{2) for q = e*'^/", 
where n is an integer greater than 2 [19]. Then for each n there are finitely 
many irreps, up to j = Each face (edge) amplitude dimj is replaced by 
dimgj, while each vertex (tetrahedron) ampUtude (6j) is replaced by {6j)q. 
In this way one obtains the Turaeev-Viro state sum invariant 

^Ty = ci,-^^E ndim,j.n(6j).W , (24) 

where dl = I]^£o(*^i™'? j)^' ^ number of tetrahedra and the {Gj)q 

symbol is appropriately normalized. 

In the case of the path integral (15), we will use the formula 

^o[r] = [lidgiU^Ilai] Wr{gi',s,,t,,) , (25) 

I f yea/ / 

where Wr is the simplical version of the spin network function for the spin 
network P. This embedded spin network is specified by labeUing a subset {I'} 
of the dual edges by the irreps {sf}, as well as by labelling the corresponding 
set of vertices {v'} with the appropriate intertwiners i^r. One can associate 
more than one irrep to a given link, in order to obtain the spin nets of the 
valence higher than four. By using (20), (21) as well as 

X,,,/^^S(^)---^SJ^)-Ec-:-eHa:5^T > ^>3 , 

(26) 
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for the edges which carry the irreps of the spin network, where C*^'^ are the 
intertwiners for n irreps, one obtains 

*o(r) = ^ n^imj/ n (6i)(w) n ^tetO/K)'^^'^)'^^''^;'^)) • (27) 

jf,H' f veV\Vr v'eVr 

The amphtude Afgt is the evaluation for the spin network consisting of the 
tetrahedral graph labelled by the jf irreps, with an additional vertex labelled 

by the intertwiner iy/, which is connected by the edges labelled by the irreps 
s to the tetrahedral graph vertices [6, 7]. This type of the spin network 
amplitudes also appears in the case of the spin foam models proposed for 
matter [21], and the spin foam amplitude (27) is a special case of a more 
general spin foam amplitude proposed for an embedded open spin network 
[21, 22]. 

As in the case of the partition function, the expression (27) is a divergent 
sum. One can then regularize it by using the SUq{2) irreps. The reason 
why this procedure works, is that the representations of SU{2) and the rep- 
resentations of SUq{2) form the same algebraical structure, i.e. the tensor 
category (see [8] for a review and references). The difference is that the 
SUq{2) category is a more general type of a tensor category, i.e. a ribbon 
tensor category, whose algebraical properties are defined by the invariance 
properties of ribbon tangles under the isotopies in or S^. The 3j and the 
6j symbols, as well as the other spin network amplitudes can be interpreted 
as morphisms in this category, i.e. the maps among the objects of the cat- 
egory and the compositions of these maps. These maps can be represented 
by the open and the closed spin network diagrams. The relations among the 
intertwiner morphisms which imply the topological invariance are preserved 
and well-defined in the quantum group case. Hence the quantum group ver- 
sion of the embedded spin network state sum (27) should be a topological 
invariant. 

Therefore a well-defined invariant for F will be given by the quantum 
group version of the spin foam amplitude (27) 

*o(r) = 51 n^™-?^/ n i^j)qi'") II Mft{jfiv'),Sl'(v'),iv',k'iv')) ■ 

jf,H' f veV\Vr v'eVr 

(28) 

The face and the vertex amplitudes are now given by the evaluations of the 
corresponding q-deformed spin networks. Nontrivial dual edge (triangle) am- 
plitudes may appear, although these can be reabsorbed into normalization 
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factors of the vertex amplitudes. In order to determine these normahzation 
factors it will be necessary to check the corresponding identities for the sim- 
plex amplitudes associated to the invariance under the motions of the graph 
in the manifold. 

4 Discussion and Conclusions 

In order for (28) to be well-defined, the spin network F should be framed, 
i.e. r should be a ribbon graph spin network. The invariant (28) could be 
interpreted as a coefficient of a state \0q) in the diffeomorphism invariant 
basis of q-deformed spin network states IF^) [1], i.e. 

|o,) = Eir.)(r<,|o,) = i:ir,)^o(r) . (29) 

The q-deformed spin network states can be defined via the q-deformed loop 
algebra [14] , and these states have the same labels as the non-deformed ones, 
but the q-deformed states have the restriction that the the maximal label j is 
given by (n — 2)/2 where q = exp{i7i/n). Since the q-dcformcd spin networks 
are naturally associated to the quantum GR with A oc 1/n [14, 1], then we 
expect that the action of the A = Hamiltonian constraint (HC) operator 
on the state (29) should give a state with the basis coefficients of order l/n^ 
or higher. Verifying this would require the knowledge of the matrix elements 
of the HC operator in the q-deformed spin network basis. 

Given that the A = quantum GR is naturally formulated in the non- 
deformed {q = 1) spin network basis [20, 1], one would like to construct a 
flat connection vacuum state in that basis. Since q — e*'^/", a natural way 
would be to consider a large n hmit of (29), i.e. 

10) = lim |0„) . (30) 

It would be interesting to study the action of the HC operator on the state 
(30), given that the matrix elements of a HC operator corresponding to the 
constraint (2) for A = are known in the non-deformed spin network basis 
[24]. 

Note that the state |0) would determine a wavefunction \E'o(^) via 

*o(^) = (^|o) = E(^|r)(r|o) = Ew^[r,^]*o(r) (3i) 

r r 

= ^Wr{gi,---,9L)M^) , (32) 
r 



9 



where g are the group elements associated to the spin network edge holonomies 
P cxp A. This expression defines an inverse loop transform. It would be 
important to understand the mathematical properties of this type of expres- 
sion, because it could serve as a definition of the vacuum wavefunction which 
was given by the heuristic expression Hies 5(e^^). Related to this, one would 
like to see what would be the result of the action of the constraint operators 
on the expression (32). 

In the case of A 7^ quantum GR, one can construct in the Tq basis a 
vacuum state 

|Ox) = Eir.)(r)i. , (33) 

whose expansion coefficients are given by the Kauffman bracket evaluation of 
the spin network F [1]. The state \0q) is different from the state \0k) because 
the invariants ^'o(r) and (r)^ are not the same. The invariant ^'o(r) is more 
complicated, since it is a sum of products of the Kaufmann brackets, which 
are the q-deformed spin network evaluations. Also, ^o(r) is by construction 
sensitive to the topology of the manifold, while (r)x is essentially an invariant 
for a spin network embedded in or R^. Since ^o(r) corresponds to A = 
quantum GR, this raises a question how to generalize the formula (33) to the 
case of arbitrary manifold topology. 

This can be done by replacing the basis coefficients of the state \0k) 
by the RT invariants for the F's (r);^^^ [5]. These invariants correspond to 
the CS theory path integrals, whose construction is based on the surgery 
representation of the manifold S and the representation theory of SUq{2). 
The RT invariant for F does not have a state-sum representation, although 
it is related to the TV invariant for F [8, 25, 26]. The TV invariant is 
given by the TV partition function for a manifold triangulation where there 
is no summation over the edge (face) irreps belonging to the embedded spin 
network. The TV and the RT invariants are related because the TV partition 
function is the square of the modulus of the RT partition function. The 
fact that the invariant ^o(r) has a very different form compared to the 
corresponding RT or TV invariants, which are naturally associated to the 
A 7^ quantum GR, is consistent with our conjecture that ^o(r) describes 
a state in the A = quantum GR. 

In order to interpret physically the vacuum states 10^) and |0) we need the 
complex connection representation and the requirement of the holomorphy 
of the corresponding wavefunction. Although a very elegant formalism, it 
is difficult to implement the reality conditions by using a scalar product. 
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We believe that this is not such a big problem, because the holomorphic 
representation seems to capture the classical reality conditions, at least in 
the case of finitely many degrees of freedom [11]. Also, the implementation of 
the reality conditions via the usual scalar product presupposes that the scalar 
products of the physical states are finite. We beheve that this assumption 
is too strong, because it is based on the interpretation of the square of the 
modulus of the wavefunction of the universe as a probability, which is a 
meaningless concept. 

Given the \0k) or the |0g) state, one is wondering is it possible to define 
an inverse loop transforms for them. On the technical side, since the algebra 
of functions on a Hopf algebra is sufficiently well understood [8], it should 
be straightforward to generalize the expression (32) to the case of functions 
on the quantum group. On the conceptual side, the resulting connection 
representation would correspond to a gauge theory based on the quantum 
group, whose properties are not well understood. 

The F — wavefunction for the A = GR has a simple form in the con- 
nection representation, i.e. \l/o(^) = One can then naturally associate 
a space of F = classical solutions to this state. However, it is clear that 
no particular F = classical solution can be associated to it, and hence we 
cannot associate this state to a flat metric vacuum because a metric flatness 
condition was never imposed, even classically. An F = classical solution 
will be a flat metric solution when A = u{e) constraint is imposed. 

Therefore if we want to construct a state which corresponds to a fiat 
metric spacetime, we would need to impose a fiat-metric condition quantum 
mechanically. In order to do this, let us consider the equation 




(34) 



In the "E" representation this equation takes a simple form 




(35) 



which has a solution 





(36) 



where A is a tensor made from the (5's such that 




(37) 
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In the "A" representation this solution will be given by 

^o{A) = jvEe'^^^'^'''^o{E) , (38) 

and let us then consider the following loop transform 

{r\Oh) ^ Jl^A W[r, A] 6{F) ^o{A) . (39) 

This suggests a fiat metric vacuum wavefunction ^o(^) = S{F)^o{A)- 
However, although it is annihilated by the constraint operators, ^'o is not 
annihilated by the A — u}{E) operator. In order to understand this, note that 
the functional integral (39) can be also interpreted as a functional integral 
over the flat connections Aq 

{r\Oh) = / VAo w[r, Ao] MAo) , (40) 

so that (39) can be interpreted as a loop transform in the partially reduced 
phase space quantization approach of Baez [9]. Since $o(^o) is annihilated 
by A — Lj{E), it would make sense to define a flat metric vacuum state as 

|o.) = Eir)(r|o.) . (41) 

r 

Given the representation (39), one could hope that the coefficients {T\Oh} 
could be defined by using the same approach as in the — 1 case, i.e. 
by using a simplical path-integral representation to define them as three- 
dimensional topologically invariant state-sums for the SU (2) quantum group 
at a root of unity. However, the main problem in this approach will be the 
non-polynomiahty of the Q{E) term, and hence one should resort to some 
sort of approximation. The simplest thing one can do is to approximate $o 
by a coherent state. For example one can take 

Q{E)= Uxj:{El-6:r , (42) 

so that $0 = e~^^^^ will be peaked around the flat metric triads e" = Sf. 
However, the expression (42) is not a diffeomorphism invariant. One can try 
to replace it by a diffeomorphism invariant analog, for example 

n{E) = d^x {det{Ei) - 1)2 . (43) 
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However, (43) is not quite a diffeomorphism invariant, since dct(£^) is a scalar 
density. Still, given a polynomial ^1{E), one can use the path- integral tech- 
niques developed for the simplical BF theories with interactions [17, 23] in 
order to obtain a state sum for a given traingulation of S. This state sum 
will not be topologically invariant if the initial D,{E) was not a diffeomor- 
phism invariant, and therefore some sort of averaging over the topologically 
equivalent triangulations would have to be performed. 

Note that this approach suggests relaxing the condition that ^o(^) should 
be annihilated by A — a) to requiring that ^'o(^) should be sufficiently peaked 
around the flat metric triads, which means using a coherent state $o- Then 
the expressions (39) and (41) also make sense in the Dirac quantization for- 
malism. 

In the case when A 7^ 0, one would need to replace the Kodama wave- 
function \]//^(yl) by ^'x(/l)$o(^) in order to construct a flat metric vac- 
uum because the de-Sitter classical vacuum solution found in [1] satisfies 
F-j — KeijkE^"' — and — 0, where is the spatial curvature, while the 
Kodama state corresponds only to F^ - — Kcij^E^"" — solutions without the 
flat metric condition. 

Matter could be included by perturbing the vacuum state via the ansatz 
^o(^) ^o(^)x[^, 0], where stands for the matter degrees of freedom, see 
[1]. In a similar manner one can also study the gravitational perturbations, 
i.e. the gravitons. Beside exploring the physical implications of the presence 
of matter, it would be also interesting to see what kind of invariants one 
would obtain by using the spin network formulation and by deflning the 
corresponding path-integrals via generalized spin foam state sums. These 
can be then compared to the generalized spin foam state sum amplitudes 
proposed for matter in [21, 22, 23]. 
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